Darboux-integration of ip = [H, f(p)] 
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A Darboux-type method of solving the nonlinear von Neu- 
mann equation ip = [H, f(p)] is developed. An explicit con- 
struction shows that self-scattering solutions are a generic 
property of such nonlinearities. An infinite-dimensional so- 
lution is constructed. 

PACS: 05.45.-a, 05.70.Ln, 03.65. Bz 



I. INTRODUCTION 



The equation 



iX = [YJ(X)}, 



(1) 



plays an important role in different branches of quantum 
and classical physics. First of all, if / satisfies /(0) = 0, 
/(l) = 1, and X = Xt = X 2 then /(X) = I | and 
the equation is simply the linear von Neumann equation 
for pure states. If X = an d Y = H is a Hamilto- 

nian, Eq. (|l]) is physically equivalent to the Schrodinger 
equation 



= H\iP) 



(2) 



for a large class of nonlinear functions /. 

Second, Eq. (@) can be written for a large class of func- 
tions f{X) as the Heisenberg equation of motion 



iX=[X,h(X)}, 



(3) 



with Hamiltonian h(X). (Assuming, for example, 
f(X) = -X n gives h{X) = X n ~ l Y + X' l ~ 2 YX + ... + 
yj" _1 .) The equations of this form are often used in 
nonlinear quantum optics. 

Third, it is known that the linear von Neumann equa- 
tion can be considered as a Lie-Poisson Hamiltonian 
system. In the case of classical "extensive" statistics, 
the Hamiltonian function is given by the average en- 
ergy Hi = Tr (Hp). Yet, there are statistical problems 
that are naturally described by nonextensive statistics 
in terms of the Tsallis q ^ 1 entropies and "average" 
energies H q = Tr (Hp q )/ Tr p q Q. One of the remark- 
able properties of q-statistics is the q-independence of 



standard geometric structures associated with equilib- 
rium thermodynamics. Extending this observation from 
equilibria to non-equilibria, one finds that Hg_is a Hamil- 
tonian function for a Lie-Poisson dynamics 



ip=[HJ(p)}, 



(4) 



where f(p) = C q (p)p q and C q (p) is a Casimir function 
satisfying C\(p) — 1. The von Neumann equation (^) 
with arbitrary nonlinearity is Lie-Poissonian with the 
Hamiltonian function Hf = Tr (if/(p)J. 
Fourth, the "Euler equations" 



iX= [Y,X% 



(5) 



appear in several contexts. The best known physical 
example (with Y kl = 5 k i/h, X mn = ie k i m Jk) is given 
by the Euler equations for a freely rotating rigid body. 
Less known and often more abstract versions of the Euler 
equations are related to the Lie-Poisson equations occur- 



ing in fluid dynamics 
abelian gauge theories 



the Nahm equations in non- 
, and the TV-wave equations for 
electromagnetic waves in nonlinear media |t|. Quite re- 
cently the equation of the form 



iX = [Y, X 3 } 



(0) 



was discovered in the context of symmetries of the Euler 
equation M. 

General equations of the form (Q) appeared in the con- 
text of nonlinear Nambu-type theories Q . Nonlinear Lie- 
Poisson density matrix equations were applied in quan- 
tum mechanics with mean-field backgrounds Jl0| and 
nonlinear quantum mechanics 0,|lJ ■ Solutions of these 
equations were used as models of non-completely-positive 
nonlinear maps which, nevertheless, satisfied physical 
conditions widely believed to be equivalent to complete 
positivity Finally the Lie-Poisson density matrix 

techniques of extending nonlinear evolutions of subsys- 
tems to entangled states proved to have applications in 
quantum information theory 

Although the literature devoted to the Euler equation 
is quite extensive ^,||,^4), analytical methods were only 
recently applied to its density matrix version p5[- p7| , p0[ . 
The constraints imposed on density matrices \(P = P-< 



1 



P > 0, Trp = 1) and Hamiltonians {W = H, H > 0, 
unboundedness) require techniques which are not based 
on standard integration via quadratures and the simi- 
lar, since the systems in question are generically infinitc- 
dimensional. The technique used in ^ [r^j20|] is an 
appropriate modification of the dressing method [ ^2|j23| 
and, in particular, of the version based on a binary Dar- 
boux transformation [|i|-^6) . 

The Darboux-type method of integration of the Euler- 
von Neumann equation 



ip=[H,p 2 } 



(7) 



introduced in ]15|] led to the discovery of the so-called self- 
scattering solutions. The process of self-scattering con- 
tinuously interpolates between two asymptotically linear 
evolutions. This very characteristic property was found 
in all the nontrivial solutions of Eq.(Q) obtained by the 
dressing method. 

A problem that remained open in all the previous 
works was how to obtain solutions of Eq.([|) with other 
values of the Tsallis parameter q. In fact, the case 
q = 2 was not very interesting from the point of view of 
nonextensive applications since the parameters involved 
in analysis of concrete physical situations were either 
close to 1 or < q < 1. Of special interest turned out 
to be the value q — 1/2 due to its plasma physics appli- 
cations. Below we present an extension of the Darboux 
technique to the nonlinear von Neumann equation for all 
functions / for which /(/?) is well defined 



II. LAX PAIRS AND THEIR DARBOUX 
CO VARIANCE 

We begin with the overdetermined linear system (Lax 
pair) 



z x {i>\ = {i>\{p-\H) 
= ~(iP\A 



(8) 
(9) 



where A, p, H are operators, A, z\ are complex numbers, 
and (ip\ is a "bra" associated with an an element of a 
Hilbert space. The operators p and H will typically play 
the roles of density matrices and Hamiltonians, respec- 
tively, but one can also think of them as just some op- 
erators without any particular quantum mechanical con- 
notations. The compatibility condition of the Lax pair 



ip=[H,A], 
[A,p}=0. 



(10) 
(11) 



It is convenient to introduce two additional Lax pairs 



z M = (p-pH)\<p), 

i\ip) = -A\<p), 
P 



(12) 
(13) 



*U<X\ = IX\{P-"H), (I 4 ) 

= ~(X\A. (15) 
v 

The method of solving ([|) is based on the following 

Theorem. Assume \ip), and (\\ are solutions of (]§), 
and <Vx|, Pi, Ax are defined by 



Pi = (l 
Ai = (l 

P = 



V — pL 

A — p 

p — V 
v 

p — V 



p 



P)P{1- 

p)a(i 



v — p 



P 



v 



P 

v — p 
P 



P 



MM 
(x\<p) ' 



Then 



z\(ipi\ 



(^(px-XH), 



(16) 
(17) 
(18) 

(19) 
(20) 

(21) 
(22) 



Proof: Substitution of Eq. (|16|) into Eq. ( |21| ) defines p\ 
by 



Pl =p+(p-v)[P,H}. 



(23) 



The proof that (23) and ( |17| ) are equivalent can be found 
in |l^] . To prove Eq. (^2|) one first notices that the oper- 
ator P given by (|19|) satisfies the nonlinear equation 



iP = -CI - P)AP - -PA(1 - P) 
p v 



(24) 



and (J2^) follows from a straightforward calculation. □ 
Lemma. If A = f(p) §, then 



ipi = [H,f( Pl )]. 



(25) 



Proof: This statement follows from the compatibility 
condition for the Lax pair (pl"|)-(p2|), if one takes into 
account 



A x = 1 



p — v 



P 



/(P)(l 



v — p 

p 



p) = f(pl)- □ 



Comment: Having one solution p one produces a new 
solution pi of the nonlinear von Neumann equation (Q). 
The procedure can be further iterated, p — > p\ — > pi — > 
. . ., in a direct analogy to the Darboux method of gen- 
erating multi-solitons |fl8| or supersymmetric quantum 
mechanics p9| . 

If p is a density matrix and v — p one can put (x\ = 
(tp\. In this case p\ is also a density matrix and spectra 
of p and pi are identical. 
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III. SELF-SCATTERING SOLUTIONS FOR A 
GENERAL NONLINEARITY 

We will now show that self-scattering solutions con- 
structed in 1 15 1 for f(p) = p 2 are a generic property of 
the general nonlinear von Neumann equation (^). 

One begins with a seed solution satisfying 



f(p) ~ a P - 



(26) 



where [A a , H] = and A a is not a multiple of the iden- 
tity. The solution satisfies 



ip=[H,f(p))=a[H,p] 



and 



p(t) 



-iaHt 



p(0)e 



iaHt 



(27) 



(28) 



Taking the Lax pair with v = p and repeating the con- 
struction from |L5| one arrives at the self-scattering solu- 
tion 



p 1 {t)=e- iam [p{Q) + {p-p)F a {t)- 1 



HI 



where 



F a (t) = { i p(0)\exp(i^--^A a t)MO)) 



(29) 



(30) 



and |<p(0)) is an initial condition for the Lax pair. 

As an example consider the harmonic oscillator Hamil- 
tonian H = u) Y^Lq n \ n ) ( n \ and 



f(p) 



2ft 



9-1 



(31) 



where q G R. To construct the operator A a we first note 
that for any q the equation 

x q - 2x q ~ 1 -x + 2 = (x"- 1 - l)(x - 2) = (32) 

has two positive solutions, x = 1 and x = 2, which will 
be used in the construction of a seed p. Choosing a = 1 
we define 



Ai = f( P ) - P = P q - Zpi- 1 - p. 



(33) 



The next task is to choose a seed p such that [p, H] ^ 
and [Ai,iZ] = 0. We do this as follows. Let |0), |1), |2), 
be the three lowest energy eigenstates of H and take 



P(0) = f(|0)(0| + |2)(2| 
7 



i(|2)(0| + |0)(2| 



(34) 



We find [p{0),H] ^ and 



A 1 =-2(|0>(0| + |2)(2| / 

-2 + |[l-(|) 1 - 9 ])|l)(l|. 



(35) 



The third eigenvalue of p(0) is 7/4 and therefore the 
density matrix is not yet normalized. Now take p = 
is/3 / (4w) . The eigenvalues of p — pH are 5/4 + iy3/4 
and 7/4 + ^-^/3/4, the latter being twice degenerated. The 
initial condition for the Lax pair is taken to be a linear 
combination of the two orthonormal eigenvectors corre- 
sponding to 7/4 + ia/3/4: 



HO)) 



V3, 



V2 



1 



V3 



|2). 



(36) 



As a general rule self-scattering solutions can occur only 
for | </5(0)) which are not eigenstates of A a . Here it means 
that q ^ 1 which is consistent with the fact that for q = 1 
the von Neumann equation is linear. 

Now we are in position to explicitly write the self- 
scattering solution for any q. We start with the seed 
solution 



p(t) = e- lHt p(0)e lHt 



and obtain 



pi(t) = e- lHt Pint (t)e 



iHt 



with 



Pint(i) = X] r { t )mn\m){n\, 



m.n—0 



rao roi r 2 
rw r n r 12 
r20 r 2 i r 2 2 



3/2 -at) at) 
-m 7/4 at) 
m m 3/2 



C(t) = 



{-3i + 73)e^* 
4-^/2(1 + e 2u} « t ) ' 
1 - i\/3 - 2e 2 ^* 
4(1 + e 2 ^*) 



(37) 
(38) 

(39) 

(40) 

(41) 
(42) 



u q = [(4/7) 1 - q - l]w/V3, uj q > for q > 1, uj q < for 
q < 1. The self-scattering asymptotics is 

^i(±oo) = 0, (43) 

C g >i(+oo) = C,<i(-oo) = -l/2, (44) 

C 9 >i(-oo) = C 9 <i(+^) - (1 - iV3)/i. (45) 

Let us note that the seed solution p(t) we have started 
with reappers in the asymtotic states 



pi — ► p, for q > 1, t — > +oo 
pi — » p, for q < 1, < —>■ — oo. 



(46) 
(47) 
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The initial condition at t = is the same for all q. 

It is interesting that, somewhat counter-intuitively, 
equations with very different q may lead to evolutions 
which are practically indistinguishable. Indeed, for q — > 
+00 we get uj q — ► oo and therefore the transition around 
t = between the asymptotic linear evolutions takes the 
less time the greater q; for q —> — oo we get uj q — ► —u>/y3. 
By the same method but with a different choice of |y(0)) 
one can generate solutions whose self-scattering takes 
place in a neighborhood of an arbitrarily chosen t = t 
(for details cf. the discussion of f(p) — p 2 given in pO|])- 

Time scales involved in self-scattering are best illus- 
trated by average position of the harmonic oscillator as 
a function of time. Fig. 1 shows (x) — Tr (xp%)/ Tr pi, 
x = (a + aJ)/-y/2 for different values of q. In Fig. 2 the 
self-scattering is explicitly seen in the contour plot of a 
Harzian ||l]] , a 3D surface representing the self-scattering 
time dependent probability density in position space as 
a function of time. 



IV. INFINITE-DIMENSIONAL EXAMPLE 

In this section we will show that the above technique 
is not limited to matrix cases. The example is, per- 
haps, rather artificial but at least clearly shows that 
a strategy leading to self-scattering solutions involving 
infinite-dimensional subspaces is possible. We are not 
interested in the trivial situation where the dynamics is 
"reducible" i.e. can be decomposed into a direct sum of 
finite-dimensional evolutions. This could be done imme- 
diately on the basis of the analysis given in the previous 
section. 

Assume that spectrum of the Hamiltonian H contains 
a discrete part {E n }™ =1 . One technical assumption we 
will need |^7j is a symmetry of the spectrum of H, at 
least when restricted to the subspace we are interested 
in. By this we mean that the spectral representation of 
H is 



H 



OO 

E 

n=i 



E n (\n,+)(n,+\ - \ n ,-){n,-\) + 



(48) 



where the dots stand for the remaining part of the spec- 
trum. 

Equation ip — [H,f(p)] is Lie-Poissonian with the 
Hamiltonian function Hj = Tr Hf(p). In nonexten- 
sive statistics one considers such generalized aver- 
ages where f(p) — p q . From the nonextensive point of 
view it is essential that it is Tr Hp q , Tr p q = 1, and not 
Tr Hp which plays the role of the energy. 

As a consequence, it is natural to treat g — f{p) nor- 
malized by Tr/(p) = 1 as a physical density matrix as- 
sociated with the solution p since then the Hamiltonian 
function has an interpretation in terms of average energy. 
In what follows we shall assume that /(— x) = f(x). 
For Tsallis-type description one can take f(x) = x q , 
q = 2n/ (2n ± f ) and the linear limit is n — > oo. 



With the above assumptions we can take the seed so- 
lution p in the form of an infinite matrix 



/ Oi 






V ! 





-oi 



«2 













-Oi 

a 3 











-«3 



\ 



(49) 



and the Hamiltonian 



H = 





Cl 





Q 








Cl 


-h 




















b-2 


Cl 














Cl 


-b 2 




















b:>, 
















C3 


-63 


V ■ 













(50) 



Any Hamiltonian whose eigenvalues are Ek = 
±y&! + |cfe| 2 can be written in this way in some basis. 
Let us note that for Ck ^ one finds [p, H ] ^ but nev- 
ertheless [g,H] — [f(p),H] = which means that the 
seed solution p is stationary. The next task is to choose 
the parameters afc and v in a way guaranteeing that the 
projector 



P 



\Xv){Xv\ 

{xAxv) 



where 



Zu{xA = (xuKp-vH) 



(51) 



(52) 



will have nonzero matrix elements between any two eigen- 
vectors of H. By construction the same will hold for 



p[l]=p + {9-v)[P,H] 



(53) 



and the nonlinear evolution will be "infinite-dimensional 
and irreducible" i.e. will involve transitions between all 
the eigenvectors of H which span the infinite-dimensional 
subspace. 

constants a, (3 satisfying bk — ua-k 
The eigenvalues of H are 
The normalization of g means 
= 1/2. It turns out that with 



Take 
and Icfe 



two 

2 



E 



± 



±a fcv /c 



! + /3 2 . 



the above choice of p the condition of "infinite dimension- 
ality and irreducibility" can be satisfied only if z v = 0. 
We look for v satisfying 



(a 2 + [3 2 )v 2 -2av+l = 0. 



(54) 



and find v± = £^jjs ■ 
The eigenvector \xv 
form 



corresponding to z v = is of the 



4 



\Xu) 



( u\w\ 

U2 W 

u k w 

V ; J 



(55) 



where 



We finally obtain 



v/2 



1 

—i 



(56) 



p[l] = 

where 

p[l] k i = a k 8 u 
G 



/ p[l] U p[l]l 2 p[l]l 3 
P[l]21 P[l]22 p[l]23 
P[l]31 P[l]32 P[l]33 



V 



1 

-1 



\ 

) 



(57) 



(58) 



a — i(3 



i 1 
1 -i 



ai 



a + i/3 



1 



i 1 



where 



2 2/3/(a„)t 



m = u k uie 



-ia[/(a*)-/(oi)]t+/9[/(o h )+/(oi)]t 



(59) 

(60) 



The solution is "infinite-dimensional and irreducible" . It 
seems that this is the first example of such a nonlinear 
dynamics one can find in the literature. The eigenvalues 
of p[l] are the same as those of p. The form of G and F k i 
shows that p[l] is again a self-scattering solution. 



V. COMMENTS 

We have shown that: (1) General von Neumann equa- 
tions with /-nonlinearities possess nontrivial solutions. 
(2) Independently of / such solutions can be obtained 
by the dressing method. (3) The self-scattering solutions 
are qualitatively similar for different nonlinearities. (4) 
The nonlinear effects are well localized in time, transient 
and asymptotically the solutions correspond to those of 
linear von Neumann equations. (5) Even large modifica- 
tions of nonlinearity may lead to small and very short- 
lived modifications of standard linear dynamics. (6) All 
nonlinearities ~ p q which are expected to be related to 
nonextesive statistics can be treated within the proposed 
formalism. (7) The formalism can be applied to genuinely 
infinite-dimensional systems. 

In the light of these findings one may wonder whether 
it is possible to experimentally distinguish between a gen- 
eral / and a linear f. Indeed, the fact that some experi- 
mental data are well fitted by a linear dynamics may only 



mean that a self-scattering has taken place in the past, 
or will take place in the future. If in addition the state 
is pure then its dynamics is given by the linear von Neu- 
mann (or Schrodingcr) equation even in case the function 
/ is highly nonlinear. It follows that not only the results 
we have reported may prove useful as a technical tool 
in many branches of classical and quantum physics, but 
they also shed new light on the negative results of exper- 
iments searching for quantum nonlinearity. 

In a separate paper we will describe a generalization 
of the above technique to a larger class of "nonabelian" 
nonlinearities |2S[] . 

Our work was supported by the Alexander von Hum- 
boldt Foundation (M.C.) and Nokia-Poland (N.V.U.). 
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FIG. 1. (x) for q = 1 (dashed), q = 2 1 / 2 (thin solid), q = n (dotted), and q = — 2 (thick solid). All evolutions for g > 1 
(g < 1) have identical asymptotic states and the same initial condition (all curves intersect at t — 0). The solution for q = 1 
satisfies the same linear equation as the asymptotic states for q / 1. 
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FIG. 2. Contour plot of a 3D surface representing probability density in position space (x\p\\x) as a function of time for 
q = 1/2, u) = 1/2, and —60 < t < 60. Continuous transition (self-scattering) between two solutions of the linear equation is 
clearly seen, pi has nonvanishing matrix elements in the subspace spanned by |2), |3), and |4). 
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